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Abstract. We present new calibrations of stellar parameters of O stars at solar metallicity taking non-LTE, wind, 
and line-blanketing effects into account. Gravities and absolute visual magnitudes are derived from results of 
recent spectroscopic analyses. Two types of effective temperature scales are derived: one from a compilation based 
on recent spectroscopic studies of a sample of massive stars - the “observational scale” - and the other from 
direct interpolations on a grid of non-LTE spherically extended line-blanketed models computed with the code 
CMFGEN (Hillier & Miller IT3351 - the “theoretical scale”. These T e fr scales are then further used together with 
the grid of models to calibrate other parameters (bolometric correction, luminosity, radius, spectroscopic mass 
and ionising fluxes) as a function of spectral type and luminosity class. Compared to the earlier calibrations of 
Vacca et al. IM the main results are: 

— The effective temperature scales of dwarfs, giants and supergiants are cooler by 2000 to 8000 K, the theoretical scale being 
slightly cooler than the observational one. The reduction is the largest for the earliest spectral types and for supergiants. 

— Bolometric corrections as a function of Tefr are reduced by 0.1 mag due to line blanketing which redistributes part of the 
UV flux in the optical range. For a given spectral type the reduction of BC is larger for early types and for supergiants. 
Typically BCs derived using the theoretical T e fr scale are 0.40 to 0.60 mag lower than that of Vacca et al. GMI, whereas 
the differences using the observational T e g scale are somewhat smaller. 

— Luminosities are reduced by 0.20 to 0.35 dex for dwarfs, by ~ 0.25 for all giants and by 0.25 to 0.35 dex for supergiants. 
The reduction is essentially the same for both T e g scales. It is independent of spectral type for giants and supergiants and 
is slightly larger for late type than for early type dwarfs. 

— Lyman continuum fluxes are reduced. Our theoretical values for the hydrogen ionising photon fluxes for dwarfs are 0.20 to 
0.80 dex lower than those of Vacca et al. (HiU, the difference being larger at late spectral types. For giants the reduction 
is of 0.25 to 0.55 dex, while for supergiants it is of 0.30 to 0.55 dex. Using the observational T e fj scale leads to smaller 
reductions at late spectral types. 

The present results should represent a significant improvement over previous calibrations, given the detailed 
treatment of non-LTE line-blanketing in the expanding atmospheres of massive stars. 
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1. Introduction 

Despite their paucity, massive stars play a crucial role in 
several fields of astrophysics: they enrich the interstellar 
medium in heavy elements; they create H II regions; they 
release huge quantities of mechanical energy through their 
winds; they explode as supernovae; they are possibly at 
the origin of Gamma-Ray Bursts; and the first massive 
stars may have reionised the early Universe at redshift 
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beyond 6. Hence, a good knowledge of their properties is 
crucial and requires the development of both evolutionary 
and atmosphere models. 

Three main ingredients have to be included in mas¬ 
sive stars atmospheres: a full non-LTE treatment since 
radiative processes are dominant over thermal processes 
(e.g. Auer & Mihalas 113721 1 : spherical expansion due to 
the stellar wind and the related velocity fields (Hamann 
119861 Hillier I1987al I1987bl Gabler et al. 119891 Kudritzki 
lu)32t Najarro et al. 119961 : and line-blanketing to take into 
account the effects of metals on the atmospheric struc¬ 
ture and emergent spectrum (Abbott & Hummer 113851 
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Schaerer & Schmutz rnroa . The two former ingredients 
were the first to be included in the models, but it is only 
recently that line-blanketing has been handled reliably. 
The main reason is the complexity of the problem which 
has to be solved when thousands of level populations from 
metals have to be computed through the resolution of 
statistical equilibrium and radiative transfer equations in 
an expanding medium. Various solutions have been devel¬ 
oped to include line-blanketing: opacity sampling method 
in the code WM-BASIC (Pauldrach et al. 1200111 , approxi¬ 
mate method to estimate line-blocking and blanketing in 
FASTWIND (Santolaya-Rey et al. 119971 Puls et al. 1200511 . 
opacity distribution functions in TLUSTY (Hubeny & 
Lanz 11995 . Lanz & Hubenv "200 21. or comoving frame cal¬ 
culations using super-levels in CMFGEN (Hillier & Miller 
urns- Each method and code has its advantages and dis¬ 
advantages: WM-BASIC makes a treatment in the ob¬ 
server’s frame uising a Sobolev plus continuum approx¬ 
imation in the solution of the rate equations and does 
not include line broadening terms, but makes a complete 
hydrodynamical calculation of the atmosphere structure; 
FASTWIND is designed for fast computations but makes 
only an approximate treatment of opacities from metals 
for which no line profile is (yet) predicted; TLUSTY makes 
a very detailed calculation of the NLTE rate equations but 
is limited to plane-parallel geometry; finally, CMFGEN 
solves the NLTE rate equations in the comoving frame 
but uses super-levels. 

The effects of line-blanketing on atmosphere models 
lead to quantitative modifications of the stellar and wind 
properties of massive stars in general and O stars in par¬ 
ticular. The most studied effect is the reduction of the 
effective temperature scale (Martins et al. 120021 Crowther 
et al. 120021 Herrero et al. 120021 Bianchi & Garcia 120021 
Repolust et al. 120041 Massey et al. 120041) . Indeed, the in¬ 
creased number of diffusions in the inner atmosphere due 
to metallic line opacities implies a heating of the deeper 
layers (backwarming effect) and consequently a higher ion¬ 
isation which shifts the relation between effective temper¬ 
ature and spectral type. The reduction can be as high as 
7000 K for extreme supergiants (Crowther et al. EES)- 
Such a change of the effective temperatures implies lower 
luminosities and lower ionising fluxes, which is crucial for 
studies involving H n regions and star forming regions. 

A good knowledge of the effective temperature scale of 
O stars is fundamental since T e g cannot be derived from 
optical photometry: the visual spectrum of O stars is in 
the Rayleight- Jeans part of the distribution and is thus al¬ 
most insensitive to T e g. Several T e g scales have been pro¬ 
posed in the past fConti ll973! Schmidt-Kalerf^HH) based 
on models without winds and metals, the most recent one 
being that of Vacca et al. csaa. As these ingredients 
are now available in models, revisions of such calibrations 
are possible. Once obtained, they can be used to calibrate 
other important parameters such as luminosities and ion¬ 
ising fluxes. 

Two main approaches lend themselves to derive a tem¬ 
peratures scale as a function of spectral type (hereafter 


ST) and luminosity class (LC): 1) the determination of av¬ 
erage T e ff’s from an observed sample of O stars, or 2) the 
determination of T e gp’s from the comparison of extended 
model atmosphere grids with the observed mean proper¬ 
ties defining the spectral types, i.e. He i/He n spectral line 
ratios for O stars. Both methods, hereafter referred to as 
the “observational scale” (1) and the “theoretical scale” 
(2) , are employed in the present work. Their advan¬ 

tages and drawbacks are roughly the following. Obviously 
the establishment of an “observational scale” relies on a 
sufficient and representative number of individual stars 
of all spectral types and luminosity classes, which have 
been analysed with state-of-the-art non-LTE line blan¬ 
keted model atmospheres with stellar winds. Several such 
studies have been pursued in the last years (Crowther et 
al. 120021 Herrero et al. 120021 Hillier et al. 120031 Bouret et 
al. 120031 Repolust et al. 120041 Markova et al. 120041 Evans 
et al. 120041 Martins et al. I2004al 120051) . However, the pa¬ 
rameter space of O stars covered by these studies based 
on non-LTE, line blanketed models with winds remains 
only partly sampled. This problem can be alleviated by 
the computation of a grid of models covering the whole 
range of parameters from which “theoretical calibrations” 
can be derived. However, as in principle the parameter 
space - accounting e.g. for M , the velocity law, microtur¬ 
bulence and others - is rather large, assumptions have to 
be made on the adopted input parameters or combina¬ 
tions thereof. Finally, these models can also be tested by 
comparison to observational results from individual stars. 
Given these advantages and limitations we have chosen 
to use both approaches to establish a new T e g scale and 
to examine their implications on the recalibration of the 
remaining stellar parameters. 

In Sect. [3 we describe the atmosphere models; we 
present the sample of observed stars from which we have 
derived estimated spectroscopic gravities and absolute vi¬ 
sual magnitudes of dwarfs, giants and supergiants in Sect. 
0 in Sect. 0 we explain our method to derive the effective 
temperature scales; in Sect.0 we show the results of the 
calibrations of other stellar parameters (bolometric cor¬ 
rection, luminosity, radius, mass, ionising flux); in Sect.0 
we discuss the effect of metallicity on the present results 
and the conclusions are gathered in Sect.0 

2. Modelling 

We have constructed a grid of models for O stars using 
the atmosphere code CMFGEN (Hillier & Miller 11998 ). 
We have tried to sample as well as possible the logg-T e ff 
plane in order to cover the whole range of spectral types 
and luminosity classes. In practice, we have relied on the 

1 Note that both scales rely on model atmospheres: one is re¬ 
lated to observations of a sample of “real” massive stars anal¬ 
ysed through spectroscopic with model atmospheres (hence the 
name “observational”), whereas the other one is only based on 
the model grid without direct link to existing stars (hence the 
name “theoretical”). 
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grid of CoStar models of Schaerer & de Koter (fTWll : mod¬ 
els A2—>E2, A3—>E3 and A4—>D4 have been recomputed 
with the same set of parameters as in the original grid. 
In addition, we have computed new models to refine this 
grid: after adopting T e g and log < 7 , we have derived log 
from the evolutionary tracks at solar metallicity of Meynet 
et al. (1994), which lead to R and then M through log < 7 ; 
the terminal velocity (uoo) was chosen to be 2.6 time the 
escape velocity (see Lamers et al. 119951 Kudritzki & Puls 
mm and the mass loss rate was computed according to 
the cooking recipe of Vink et al. (12001 1) . A total of 38 mod¬ 
els with 3.0 < log g < 4.25 and 25000 < T e g < 48000 K 
was finally computed. 

The main characteristics of CMFGEN are largely de¬ 
scribed by Hillier & Miller 119981) . We give here a short 
overview: 

- Non-LTE: all statistical equilibrium equations are 
solved individually to provide the detailed level popula¬ 
tions. The radiative transfer equation is also solved and 
leads to the radiation field. 

- Line-blanketing: line-blanketing is included directly 
by means of a super-level approach. In addition to H and 
He, C, N, O, Si, S and Fe level populations have been 
computed in our models. This amounts to ~ 2000 lev¬ 
els and ~ 800 super-levels. Solar abundances according to 
Grevesse & Sauval (T1W8D have been adopted. Numerous 
tests have been undertaken to confirm the accuracy of 
the super-level approach. While the strength of individ¬ 
ual lines can be affected by using too few super-levels, 
or by poor super-level assignments, the method provides 
consistent non-LTE blanketed structures. Moreover, the 
super-level assignments are easily changed, and are often 
simply used to group levels of the same LS configuration. 
Note that the adopted super-level approach does not di¬ 
rectly influence line blanketing: its influence is indirect, 
and arises through its effect on the non-LTE level popu¬ 
lations obtained from the solution of the statistical equi¬ 
librium equations 

- Temperature: the temperature structure is derived 
from the radiative equilibrium equation. 

- Hydrodynamical structure: CMFGEN does not com¬ 
pute self-consistently the hydrodynamical structure which 
has to be given as input and which is held fixed. We have 
adopted the velocity structure resulting from the smooth 
connection by smooth connection, we mean that both 
the velocity v(r) and its derivative dv/dr are imposed 
to be continuous at the connecting point, of a pseudo¬ 
hydrostatic photospheric velocity structure with the /3-law 
(v = Vqq (1 — |r)^) representing the wind (with a value of 
/3 chosen to be 0.8). The photospheric structure was taken 
from the TLUSTY grid of models OSTAR2002 (Lanz & 
Hubenv f20()2jl . Fig.^shows the initial TLUSTY structure 
(dashed line) together with the CMFGEN input structure 
(solid line) in a typical case. One sees that the transi¬ 
tion between the TLUSTY hydrostatic structure and the 
CMFGEN “wind” structure corresponds to a velocity of 
the order of a few km s -1 (or typically 10 to 30 % of 


the sound speed) which ensures that the important effects 
of velocity fields and spherical expansion are taken into 
account even in the photosphere (e.g. Kudritzki nrn . 
Note also that the density at the connecting point scales 
with the mass loss rate, so that the adopted wind struc¬ 
ture has most effects in supergiant models. When the pa¬ 
rameters T e ff-log <7 did not match any of the TLUSTY 
models, we have simply linearly interpolated the photo¬ 
spheric velocity structure from the closest neighbours in 
the OSTAR2002 grid (see also Martins et al. mm . For ex¬ 
ample for a model with temperature 35000 < T e g < 40000 
K, we estimated the density and optical depth from the 
TLUSTY models with T e g = 35000 and 40000 K accord¬ 
ing to x = 4000(^-35000 * (*40 — ir35) + a:35 where xi is 
the density or optical depth at temperature i in kK. This 
method, although approximate, was tested by Bouret et 
al. (120030 and was found to give the same structure as a 
true computation with TLUSTY. 

- Microturbulent velocity: a microturbulent velocity 
can be included in the models. In the computation of the 
atmospheric structure (level populations, radiation field), 
we have chosen a value of 20 km s _1 while in the computa¬ 
tion of the formal solution of the radiative transfer leading 
to the detailed emergent spectrum, we have adopted Vturb 
= 5 km s -1 . The quite high value used in the atmospheric 
structure is required to avoid extremely long computation 
time. Test models have revealed that the He I A4471 and 
He II A4542 line profiles were basically unaffected by a 
decrease of r> tU rb from 20 km s -1 to 10 km s -1 in the at¬ 
mosphere structure calculation (see Martins et al. I2004al 
Martins et al. mm- Note also that microturbulent veloc¬ 
ities of 10 to 20 km s ^ 1 are often found when modelling 
supergiants (e.g. hillier et al. mm- 

3. Observed sample 

We have compiled the results from optical spectroscopic 
analysis of Galactic O stars using non-LTE spherically ex¬ 
panding models including line-blanketing. In practice, we 
have used the results of the studies of Herrero et al. ma 
and Repolust et al. fgnna based on models computed with 
FASTWIND. We have also used our own results from the 
analysis of Galactic O dwarfs with CMFGEN (Martins et 
a imm so that the number of O stars in this observed 
sample amounts to 45 objects. We have not included the 
results of Markova et al. since most of the stellar 

parameters where derived from calibrations and not from 
analysis with atmosphere models. We also excluded ob¬ 
jects analysed by Bianchi & Garcia (120021) and Garcia & 
Bianchi cssa since their studies are based on pure UV 
analysis and yield quite discrepant results to optical and 
combined UV-optical analysis for reasons we discuss now. 

The above mentioned pure UV analysis (using WM- 
BASIC) relied on two main T e g indicators, O V A1371 
and P V AA1118,1128 , lines which strongly depend on 
the wind properties (in particular clumping, see Crowther 
et al. nia Bouret et al. mm- Moreover, WM-BASIC 
does not allow a treatment of Stark broadening which 
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Fig. 1. Typical hydrodynamical structure: density as a 
function of Rosseland optical depth for a model with T e g 
= 40000 K, logg = 3.50 and M = 1.410 -5 M 0 /yr. The red 
dashed line is the initial TLUSTY structure on which a (3 
velocity law is connected to give the final CMFGEN in¬ 
put structure (black solid line). The dotted line shows the 
CMFGEN velocity structure as a function of Rosseland 
optical depth. The TLUSTY and CMFGEN structures are 
similar in the very inner atmosphere and start to depart 
from each other around log TR OSSe i an d ~ 0 which corre¬ 
sponds to a velocity of a few km s _1 . 

is crucial to correctly predict the width of gravity sen¬ 
sitive lines which are moreover essentially absent in the 
UV range. Note that the effects of clumping should not 
directly modify the present results, since we use photo- 
spheric lines formed well below the radius where clumping 
appears. In spite of this, we can not discard the fact that 
UV T e ft diagnostics sometimes lead to lower values than 
optical lines, as is the case of the WM-BASIC studies of 
Bianclii & Garcia 120021) and Garcia & Bianchi n- 
However, other recent analysis show that consistent re¬ 
sults with good fits can be obtained for photospheric UV 
and optical lines (Bouret et al. |20()(>1 Hillier et al. 2'(')fi3 
Martins et al. vm. Given this, and since our present 
approach is based on T e g- sensitive optical He lines, we 
do not include the results of spectroscopic analysis based 
only on UV non-photospheric lines in our observational 
sample. 

4. Gravities and effective temperatures 

As gravity is the main parameter determining the lumi¬ 
nosity class a calibration of this quantity is needed to allow 
the calibration of effective temperature for different LC. 


Once established, such a log g-ST relation for all LC can 
be used in the model atmosphere grid to produce new re¬ 
lations between effective temperature and spectral type, 
using gravity to discriminate between luminosity classes. 
These steps, as well as the determination of the observa¬ 
tional T e s scale are now discussed. 

All the resulting calibrated stellar parameters are sum¬ 
marised in Tables n t° El 

4.1. log g-ST relation 

The observational sample including recent spectroscopic 
analysis of individual stars has been used to derived an 
empirical calibration log g-ST. Note that the gravities 
used here are corrected for the effect of rotation, except for 
the dwarfs stars studied by Martins et al. HIB . However 
in the latter study, the rotational velocities are not larger 
than 130 km s -1 so that the difference between derived 
logg and “true” logg is negligible (see e.g. Repolust et 
al. I2004jl . The results are shown in Fig. |2 where dwarfs 
are represented by triangles, giants by squares and super¬ 
giants by circles. We note that there is a significant scatter 
in the empirical points leading to standard deviations of 
0.15, 0.07 and 0.12 dex respectively for dwarfs, giants and 
supergiants. The linear fits are to the observational results 
for dwarfs, giants and supergiants are: 

f 3.924- 0.001 x ST (V) 

logg = l 3.887 - 0.040 x ST (III) (1) 

[ 3.979- 0.083 x ST (I) 

where ST is the spectral type. 

Vacca et al. CEUi established a relation between spec¬ 
troscopic gravities and spectral type from results of spec¬ 
troscopic analysis through plane-parallel pure H He non- 
LTE models. We have proceeded similarly except that we 
have relied on atmosphere models including also winds 
and line-blanketing. Figure^ (upper panel) shows the dif¬ 
ferences between the two calibrations for various luminos¬ 
ity classes. We see that they are very close, the difference 
being less than 0.1 dex, which is the typical error gener¬ 
ally quoted for the determination of spectroscopic grav¬ 
ities from the fit of Balmer lines. Herrero et al. (120001) 
showed that the inclusion of winds in atmosphere mod¬ 
els lead to a systematic increase of logg by 0.1 to 0.25 
dex. The present results show that on average gravities de¬ 
rived with plane parallel H He models are similar to those 
derived with line blanketed spherically extended models. 
Note that the results of Herrero et al. (120001) are for su¬ 
pergiants with strong winds (M of the order of 10 -5 "~ 6 
M 0 /yr) for which an effect of the wind on the log g sen¬ 
sitive lines is not unexpected. For stars with such high 
density winds, the lowering of T e ff imposed by the inclu¬ 
sion of line-blanketing may compensate for this wind effect 
so that, fortuitously, the relation logg - ST is not signifi¬ 
cantly changed. 
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Fig. 2. Lower panel : gravity as a function of spectral type. 
Symbols corresponds to results derived from spectroscopic 
analysis of observed Galactic stars (sources: Repolust et 
al. [27TO1 Herrero et al. mm and Martins et al. m . 
Triangles are dwarfs, squares are giants and circles are 
supergiants. Lines are the least square fits of the obser¬ 
vational data (solid: dwarfs; dashed: giants; dot-dashed: 
supergiants). Upper panel: difference between the calibra¬ 
tion of Vacca et al. mm and the present ones. 

4.2. T e ffscales 

As already discussed, the derivation of new effective tem¬ 
perature scales was made following two different meth¬ 
ods: First, the relation between spectral type and T e g 
was derived directly from our grid of models, using the 
above log g— ST relation to discriminate between luminos¬ 
ity classes (“theoretical scale”). Second, we used the sam¬ 
ple of stars analysed with quantitative spectroscopy to es¬ 
timate the mean relation T e g-ST (“observational scale”). 
The resulting calibrated stellar parameters using the the¬ 
oretical T e ff scale are summarised in Tables U to El those 
based on the observational T e ff scale in Tables IS to 0 

4.2.1. T e g scale from models (“Theoretical scale”) 

Here, we show how we derived new T e g scales from our 
grid of models. The advantage of using models is that 
we can sample the whole spectral and luminosity range. 
However, the results depends on the grid of model and the 
parameters adopted to compute these models. 

The basic idea is to build a function T e g (ST,log g) from 
which temperatures from dwarfs, giants and supergiants 
— characterised by a given relation logg-ST — were de¬ 
rived. For this, we have made 2D interpolations using the 


NAG routines eOlsef and eOlsff. The basic principle of 
such routines is to construct a surface T e ff(logg,ST) regu¬ 
larly sampled in log g-T e g from the grid model points (not 
regularly sampled in the space parameter) and then to 
determine the effective temperatures from this surface for 
a given log g-ST relation. The construction of the surface 
involves two parameters representing basically the num¬ 
ber of neighbouring input points from which we want the 
interpolation to be made for a given (log g,T e fi) of the 
regular grid. We have tested the effects of a change of 
the values of these parameters on our results and have 
found that for a reasonable number of neighbouring points 
(> 15) there was essentially no modification. Note that 
we have not used directly the spectral type as a param¬ 
eter, but we rely on the ratio of the equivalent widths 
log(EW (H el X4471) / EW (He//A4542)) = log W'), which 
is directly related to ST through the classification scheme 
of Mathys (119881) . According to this scheme, a given spec¬ 
tral type corresponds to a range of values for W' so that 
all stars having W' in this range will have the same spec¬ 
tral type, introducing thus a dispersion in any calibration 
as a function of ST. This interpolation process leads to the 
T e ff-ST relations displayed in Fig. El The 1 a dispersion 
on the logg-ST relation translates to an error of ±500A' 
for the late spectral types and up to ±1000A' for early 
spectral types, with little differences between luminosity 
classes. Note that between spectral types 09.5 and 09.7 
there is only 0.2 unit in terms of spectral type, but the 
difference in terms of log W' is 0.3 dex, which is as large 
as between spectral types 07.5 and 09. Hence, we have 
excluded the 09.7 spectral type in our calibrations since 
the interpolations (made on loglW and not on spectral 
type) lead to artificially low effective temperatures. 

The differences between our temperature scales and 
the calibrations of Vacca et al. (119901) are plotted in the 
upper panel of Fig. El Typically they reach from 2000 and 
8000 K, the largest differences being found for supergiants. 

We note also that our effective temperatures for the 
latest spectral types (09.5) agree reasonably well with 
the effective temperatures of the earliest B stars. Indeed, 
Schmidt-Kaler 119821) predict T e g = 30000 K (29000, 
26000) for dwarfs (giants, supergiants) of spectral type 
BO, while we derive T eff = 30500 K (30200, 28400) for 
09.5 dwarfs (giants, supergiants). 

4.2.2. “Observational” T e g scale from a sample of 
Galactic O stars 

Instead of relying on a grid of models, one can simply de¬ 
rive mean T e g scales from the results of detailed spectro¬ 
scopic analysis of O stars (see also Repolust et al. |20'0"I| ). 
The advantage is that the parameters of the best models 
are hopefully the true stellar and wind parameters of the 
stars analysed. But the problem is that the number of ob¬ 
jects is limited and that all spectral types are not available 
(this is especially true for giants). In Figs. El Eland El we 
show the results of a linear regression to the data points 
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Fig. 3. Lower panel: “theoretical” effective temperature 
scale for dwarfs (solid line), giants (dashed line) and super¬ 
giants (dot-dashed line) as derived from our model grid. 
Upper panel: difference between the T e g scales of Vacca et 
al. (1199611 and the present results. 

(dashed line) in comparison to the T e g scale derived from 
the grid of models (solid line). 

The linear fits to the observed data are: 

( 50838 - 1995 x ST (V) 

T eff = l 50882 - 2115 x ST (III) (2) 

[ 47666 - 1881 x ST (I) 

and the dispersions are 1021, 529 and 1446 K respectively. 

The differences with the relations of Vacca et al. (119961) 
are also shown in the upper panels of the figures and are 
slightly smaller (2000 K) than for the theoretical scale, 
mainly for late spectral types. 

4.2.3. Theoretical vs observational T e g scale 

At this point, it is worth discussing the differences between 
the “theoretical” and the “observational” relation. First of 
all, given the uncertainties and the natural dispersion of 
the observational results (~ 2000-3000 K) the agreement 
between our theoretical and observational scales is good 
for dwarfs and supergiants, although our effective temper¬ 
atures seem to be slightly underestimated (by ~ 1000-2000 
K) for late spectral types in our models. For giants, our 
theoretical T e g scale seems to be a little too cool (by ~ 
1500 K). Note, however, that the number of giants studied 
so far is relatively small (7). 

The following assumptions could be at the origin of 
the lower effective temperatures in our models compared 
to the results from the analysis of individual O stars: 1) 



2 4 6 8 10 

Spectral Type 


Fig. 4. Effective temperature scale for dwarfs as derived 
from our model grid (“theoretical” scale, solid line) com¬ 
pared to “observational” results based on non-LTE line 
blanketed models including winds (Herrero et al. . 2002 
Repolust et al. 120041 and Martins et al. 120051 ). The dashed 
line is the result of linear regression to the data points. 
The typical error bars for the observational results are 
±0.5 unit for spectral types and are taken from detailed 
spectroscopic analysis for T e g. The upper panel shows the 
difference between the relation of Vacca et al. CM and 
our T e ff scales (solid line: relation from model grid; dashed 
line: relation from spectroscopic analysis). 

incorrect metallicity/abundances, 2) incorrect atmosphere 
and wind parameters (microturbulent velocities, mass loss 
rates, slope of the wind velocity law), or 3) incorrect grav¬ 
ities. We now discuss these one by one. 

1) Metallicity/abundances: First, it is well estab¬ 
lished that line-blanketing leads to lower T e g so that a too 
strong decrease could be due to a too high metal content. 
However, all stars used for the comparison are Galactic 
stars which should have solar composition. Of course there 
is a metallicity gradient in the Galaxy, but the metal con¬ 
tent remains within approximately ± 2 times the solar 
content (e.g. Giveon et al. I2002a[) . Since we have shown 
that increasing the effective temperature by 1000 to 2000 
K required a decrease of the metal content by a factor of 
8 (paper I and Martins I2004bl) . it is not likely that the 
difference between our new T e g scales and observational 
results can be explained by different metal contents. 

2) Atmosphere and wind parameters: Another 
possibility is that some atmospheric and wind parameters 
used for our modelling of O stars atmosphere are not ex¬ 
actly the same as the parameters of real O stars. In paper 
I, we have tested a number of such parameters (microtur- 
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Fig. 5. Same as Fig. 0for giants. 



Spectral Type 

Fig. 6. Same as Fig. U| for supergiants. 

bulent velocity, mass loss rates, slope of the velocity field in 
the wind -/? parameter-) and shown that increasing z> tur b, 
decreasing (3 or decreasing M could lead to later spectral 
types for a given T e g (or equivalently to higher T e g for 
a given spectral type). Since (3 = 0.8 is already small in 
our models, it is not likely that this parameter can ex¬ 
plain the difference. However, mass loss can explain the 
discrepancy since most of the late type observed dwarfs 
in Fig El are from Martins et al. mm and have M much 


lower than the predictions of the Vink et al. II2001D cook¬ 
ing recipe used to assign a mass loss rates in the present 
models. Indeed, Martins et al. found that dwarfs 

with luminosities below ~ 10 52 L 0 have lower mass rates 
than the results of hydrodynamical simulations. But Fig. 
^ shows that the discrepancy between our derived effec¬ 
tive temperatures and the observed ones starts for stars 
with spectral types later than 06, and it turns out that 
according to our new calibrations an 06V star has a lumi¬ 
nosity of ~ 10 5 ' 3 L 0 . We have also shown in paper I that 
increasing M by a factor 2 lead to a decrease of logW' 
by ~ 0.1 dex. Hence, it is likely that mass loss is at the 
origin of the underestimation of T e Q in our models of late 
type dwarfs compared to results of spectroscopic analysis. 
Concerning supergiants, the situation is a little more com¬ 
plicated since according to Tabled all supergiants have lu¬ 
minosities higher than 10 5 ' 2 L 0 and should therefore not 
have reduced winds. Hence, the reason for the lower T e g in 
our models is probably not rooted in the mass loss rates. 

Mokiem et al. rcnna have also shown, from the com¬ 
putation of theoretical spectra, that an increase of micro- 
turbulent velocity from 5 to 20 km s _1 translates to a shift 
towards later spectral types by nearly half a sub-type since 
both He I and He II lines are stronger but the former are 
more affected (cf. paper I), although the effect is usually 
more important on singlet than on triplet He I lines (see 
also Smith & Howarth ll998l Villamariz & Herrero 2000 )). 
This can partly explain the observed behaviour for super¬ 
giants for which microturbulent velocities of 20 km s -1 
are possible. Test computation in our models reveal that 
increasing i> tur b by such an amount corresponds to a shift 
in loglW by ~ 0.15, which is of the order of the width of 
the range of W' values within a spectral type. Note also 
that in our grid of models, we have used n tU rb = 20 km 
s” 1 to compute the atmospheric structure. This may be a 
little to high for dwarfs, but as already discussed (Sect.0 
test models reveal that this does not lead to an overesti¬ 
mation of the line-blanketing effect (only some particular 
lines are modified). 

3) Gravity: Another possibility to explain the differ¬ 
ent T e ff - scales is that gravities from spectroscopic analy¬ 
sis may be underestimated. This would lead to lower gravi¬ 
ties for a given luminosity class, and consequently to lower 
effective temperatures. However, if we recompute the T e g 
scale with gravities increased by the standard deviations 
obtained from least square fits of observed data, the corre¬ 
sponding increase in T e g we find (~ 1000 K, see Sect. ro 
is not sufficient to explain the difference with observed 
effective temperatures. This is especially true for giants 
for which the theoretical effective temperatures are sys¬ 
tematically underestimated by ~ 1500 K. Note however 
that the small number of giants studied so far may hide 
the natural spread in the T e g - spectral type relation we 
see for dwarfs and supergiants. This may thus artificially 
increase the difference with our theoretical results. 

We also want to mention that we have not included the 
results of Bianchi & Garcia J2303 and Garcia & Bianchi 
(12001 in our comparisons since they are based on UV 
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Fig. 7. Comparison between H and He line profiles of the 
same model (A2, see paper I) computed with the Isa-Wind 
photosplreric structure (dashed line) and the TLUSTY 
structure (solid line). See text for discussion. 

analysis, whereas both the present modelling and the ob¬ 
servational results of Repolust et al. Herrero et 

al. 01717211 and Martins et al. (1200511 rely on optical di¬ 
agnostics. Bianchi & Garcia rcnna and Garcia & Bianchi 
rcnna found effective temperature much lower (by several 
thousands K) compared to any optical study. This may 
be partly due to the dependence of the UV diagnostics 
on several parameters other than T e g , especially clump¬ 
ing which may significantly alter the ionisation structure 
in the wind (but not the photosphere) and thus the 
strength of UV lines and which is not taken into account 
in their analysis. 

We note that the present T e g scale for O dwarfs is 
slightly cooler (~ 500 K) than that of paper I. This is ex¬ 
plained first by the lower gravities used here for dwarfs 
(values of 4.00-4.05 were adopted for logg in paper I) 
and second by the different pseudo-hydrostatic structure 
used: ISA-Wind structures were computed in paper I while 
TLUSTY structures are used here. The difference is illus¬ 
trated in Fig. 0 where we see the effect of the photospheric 
structure on the He classification lines in a model with 
T e g = 33300 K. TLUSTY structures, which are more re¬ 
alistic due to a better treatment of line-blanketing, have 
higher temperature and thus a higher thermal pressure 
which implies a higher density (through the hydrostatic 
equilibrium equation). This leads to slightly higher ion¬ 
isations and then to earlier spectral types. The present 
results should represent an improvement over paper I. 

The general conclusion of the above discussion is that 
our theoretical T e g scale should be taken as indicative and 


have an uncertainty of ± 1000 to 2000 K due to both a 
natural dispersion and uncertainties inherent to our global 
method. Hence in the following, we use both these theo¬ 
retical T e s scales AND the relation derived from detailed 
modelling of individual stars (“observational” relations) 
to calibrate the other stellar parameters. 

5. Other stellar parameters 

We now derive calibrations for the remaining principal 
stellar parameters based on the two temperature scales 
determined above. The results are summarised in Tables 

mto@ 

5.1. Magnitudes, bolometric corrections, and 
luminosities 

The calibration of luminosity requires two ingredients: 
first the absolute visual magnitude and second the bolo¬ 
metric corrections. To calibrate My as a function of spec¬ 
tral type, we can either rely entirely on models and pro¬ 
ceed exactly as for the effective temperature, i.e. adopting 
the logg-ST relations, or we can derive empirical calibra¬ 
tions of My from results of analysis of individual stars, 
as we have done for logg. Here, we have chosen the latter 
method. The advantage is that the derived My-ST cal¬ 
ibration depends directly on observations and relies less 
on atmosphere models. It may however suffer from the 
often poorly know distances of the analysed stars, render¬ 
ing more uncertain the absolute magnitudes. The former 
method (log g-ST + My from models) relies on better My 
but suffers from the uncertainty in the logg-ST calibra¬ 
tion. As the final goal is to produce calibrations as close 
as possible to observations, we do not retain this method. 

The result of the empirical calibration My-ST is shown 
in Fig. |S1 where the symbols and observational data are 
the same as in Fig. |3 The standard deviation is 0.40 for 
dwarfs, 0.26 for giants and 0.45 for supergiants. The linear 
fits to the observational data are: 

f 0.291 x ST - 6.667 (V) 

M v = l 0.147 x ST - 6.570 (III) (3) 

[ 0.011 x ST - 6.388 (I) 

where ST is the spectral type. One may argue that due 
to the fact that massive stars evolve almost horizontally 
in the HR diagram after the main sequence, there could 
be a wide range of luminosities (and My) for a given T e g 
(or equivalently ST) for supergiants so that a calibration 
My - ST for luminosity class I stars may be questionable. 
However, the observed dispersion observed in the present 
study and in the Vacca et al. (M study is not signif¬ 
icantly different between dwarfs, giants and supergiants. 
Also, the observed samples show a clear separation in My 
between the different luminosity classes which is even more 
highlighted by the fact that we consider only 3 luminos¬ 
ity classes (and not 5 as usual). Hence, we think that a 
calibration My - ST even for supergiants makes sense. 
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Fig. 8. Same as Fig.^for absolute magnitudes. The stan¬ 
dard deviation from the mean relation is 0.4 mag for 
dwarfs, 0.26 mag for giants and 0.45 mag for supergiants. 

Absolute visual magnitudes have also been derived 
from observational results by a number of authors. The 
upper panel of Fig. |H1 shows the difference between the 
calibrations of Vacca et al. and the present ones. 

They do not differ by more than 0.4 dex which is also 
the typical dispersion of our relations. Hence, we can not 
quantitatively compare both types of relations: agreement 
or small differences (< 0.4 dex) can not be distinguished. 

Recently, Schroder et al. JW)31l have tested the My- 
ST relations of Schmidt-Kaler mm, Howarth & Prinja 
imi and Vacca et al. by comparing the parallaxes 

derived from such relations to the parallaxes measured by 
Hipparcos for a number of Galactic O stars. They con¬ 
cluded that the three calibrations are consistent with the 
Hipparcos measurement. As our relations are not strongly 
different from the relations of Vacca et al. iTTT^ tested 
by Schroder et al. mm . they are certainly reasonable. 

To obtain the luminosities from absolute magnitudes, 
one needs to know the bolometric correction. Vacca et al. 
ITSSEI have shown that bolometric corrections were essen¬ 
tially independent on the luminosity class and depended 
mainly on T e ff. We have computed the BC of our models 
from luminosities and synthetic My, yielding the following 
linear dependence on T e g: 

BC = 27.58 - 6.80 x log T eS (4) 

with a standard deviation of 0.05. Our relation is shifted 
by 0.1 magnitude compared to that of Vacca et al. 
as seen on Fig. [9] This is explained by the effect of line¬ 
blanketing which redistributes the flux blocked at short 
wavelengths to longer wavelengths, including the optical. 



Fig. 9. Bolometric correction as a function of effective 
temperature from our set of CMFGEN models. The solid 
line is the linear regression (Eq. and the dotted line is 
the relation of Vacca et al. and is shifted by ~ 0.1 

mag compared to our relation. 

This translates to an increase of the V magnitude by ~ 
0.1 mag, and consequently to a reduction of the bolomet¬ 
ric correction by the same amount (for a model of given 
luminosity). 

Figure fTiH shows our calibration of bolometric correc¬ 
tions as a function of spectral type obtained with our the¬ 
oretical T e ff scale together with the relations of Vacca et 
al. H1996|l . The difference between both types of relations 
is displayed in the upper panel. The reduction of BC goes 
from ~ 0.4 to ~ 0.5 magnitude for dwarfs, from ~ 0.35 
to ~ 0.6 mag for giants and from ~ 0.4 to ~ 0.65 mag 
for supergiants. Figure shows the bolometric correc¬ 
tions estimated with both the “theoretical” and “observa¬ 
tional” T e ff scales. Both types of BCs differ by less than 
0.2 magnitudes (which is lower than the difference with 
the results of Vacca et al. 1199611 . except for the latest su¬ 
pergiants. The lower BCs at late spectral types for the 
theoretical relation simply comes from the lower effective 
temperatures given by the theoretical T e g scales. 

With the absolute magnitudes and bolometric correc¬ 
tions, the luminosities can be estimated according to 

log -f- = -0.4 x (My +BC- Mg 1 ) (5) 

Lq 

where Mg 1 is the solar bolometric luminosity taken to be 
4.75 (Allen |T57B|. Again, both theoretical and observa¬ 
tional T e ff scales can be used. To estimate the error on 
the luminosities, we have recomputed them adding ±cr on 
logg and My in the logg-ST and My-ST relations for 
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Fig. 10. Bolometric correction as a function of spectral 
type from our set of CMFGEN models (solid lines) and 
from Vacca et al. C225I dotted lines). BCs are usually 
reduced by ~ 0.2 to 0.5 magnitudes for a given spectral 
type, the reduction being the largest for early supergiants. 



Fig. 11. Same as Fig. 1101 but with the bolometric cor¬ 
rections derived using the “theoretical” (bold lines) and 
“observational” (thin lines) T e g scales. The upper panel 
shows the differences between “theoretical” and “observa¬ 
tional” BCs. 


a given T e g scale. This means that these errors do not 
take into account the uncertainty on the effective tem¬ 
peratures. This results in errors of ±0.15,0.10 and 0.17 
in log jf- for dwarfs, giants and supergiants respectively. 
F igure IT?1 shows the resulting average positions of dwarfs, 
giants and supergiants in the HR diagram. For a given 
spectral type, the luminosities are lower by 0.25-0.35 for 
dwarfs, by ~ 0.25 dex for giants and by 0.25 to 0.35 dex for 
supergiants compared to the values presented by Vacca et 
al. (1199611 . Given the small differences in BC for the two 
T e g scales it is clear that L depends little on this choice. 
Quantitatively, using the theoretical or observational T e g 
scales results in modifications of the luminosity by less 
than 0.1 dex for all LG. 

A comment concerning the off-set of the dwarfs posi¬ 
tion compared to the Zero Age Main Sequence seen in Fig. 
H2lis worthwhile. As we are determining the average prop¬ 
erties of O stars, it is not surprising that the O dwarfs rela¬ 
tion does not fall on the Zero Age Main Sequence (ZAMS). 
However, a precise explanation of the offset, indicating av¬ 
erage ages of ~ 1 to 5 Myr for early to late O dwarfs, is still 
lacking. A possible explanation could be that very young 
stars still embedded in their parental cloud may not be 
observable. As later types have weaker winds less efficient 
to blow this cloud, this could possibly even explain why 
later O dwarfs have older average ages. However, although 
few, some O dwarfs are found very close to the ZAMS 
(Rauw 120031 Martins et al. 120051) . Systematic studies of 
larger samples of the youngest O stars will be needed to 


settle this question. Finally, the existence of a theoreti¬ 
cal ZAMS is not fully established (Bernasconi & Maeder 
11996 ) and the precise shape/location of the ZAMS and 
isochrones is also dependent on parameters like the rota¬ 
tion rate (Meynet & Maeder I2000JI . 

5.2. Mass and radius 

With the effective temperature scale and the luminosities 
as a function of spectral type, it is easy to derive the re¬ 
lation radius - spectral type from 


R = 


4ttct ijT e 4 ff ’ 


( 6 ) 


where an is the Stefan-Boltzmann constant. The standard 
deviation on R is 10 to 20 %. 

Gravity can then be used to derive the spectroscopic 
mass according to 



(7) 


with G the gravitational constant. Masses must be taken 
with care since their uncertainty is as high as 35 to 50 %. 
For all results see Tables 12 to El 
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Fig. 12. HR diagram showing the typical position of 
dwarfs (triangles), giants (squares) and supergiants (cir¬ 
cles) from the present “theoretical” calibration of T e g and 
log . Isochrones for 0, 1, 3 and 5 Myrs are shown to¬ 
gether with evolutionary tracks for different masses (taken 
from Meynet et al. GHH1 . The crosses and dotted lines are 
the positions of dwarfs, giants and supergiants following 
the Vacca et al. (119961) calibrations. 

5.3. Ionising fluxes 
5.3.1. Calibrations 

An important quantity related to massive stars is the num¬ 
ber of ionising photons they emit due to their high effective 
temperature and luminosity. Ionising fluxes are defined as 
follows 

[ Xi *XFx„ 

qi= —— d\ (8) 

Jo he 

where F\ is the flux expressed in erg/s/cm 2 /A and A,; is 
the limiting wavelength below which we estimate the ion¬ 
ising flux. In practice, we usually define go, Qi and g 2 as 
the fluxes able to ionise Hydrogen (Ao = 912 A), He I (Ai 
= 504 A) and He II (A 2 = 228 A). The ionising fluxes 
are highly sensitive to line-blanketing effects since metals 
have most of their bound-free transitions in the Lyman 
continuum. Hence, the presence of such bound-free opac¬ 
ities in metals modifies the spectral energy distribution 
(SED) and consequently the ionising fluxes: flux at high 
energy is blocked and must be re-emitted at longer wave¬ 
length. Moreover, the numerous metallic lines in this range 
create a veil which can also significantly alter the shape 
of the SED. 

Several studies have already provided new ionising 
fluxes as a function of effective temperature (Smith et al. 


Mokiem et al. 12004 1. An in depth study of ionising 
fluxes as a function of spectral type and metallicity (and 
thus line-blanketing) was carried out by Kudritzki 1 2002 1 
who found that the number of Lyman continuum photons 
was not a strong function of metallicity. This may be sur¬ 
prising at first glance since, as we have already mentioned 
in Sect. EU flux redistribution takes place at all wave¬ 
lengths, including A > 912 A. However, energy flux and 
photon flux are not equivalent. It is possible for the flux in 
the Lyman continuum to go down while the photon flux 
stays almost constant. In this case, most of the photons re¬ 
side close to the Lyman edge while energy is redistributed 
at all wavelengths, including A > 912 A. 

In the studies mentioned above, ionising fluxes where 
calibrated as a function of effective temperature but not 
as a function of spectral type. Due to the cooler new effec¬ 
tive temperature scales, such calibrations must be signif¬ 
icantly different from previous relations based on unblan¬ 
keted models. For astrophysical applications, the interest¬ 
ing quantity is in fact Qi, the total number of ionising 
photons emitted per unit time related to qi by 

Qi = q t x 4-7T R 2 , (9) 

with R the stellar radius. With our calibrations T e ff~ST 
and R-ST, we are now able to derive relations between 
spectral types and the total ionising fluxes Qi. However, 
we have not considered the He II ionising fluxes since they 
may be strongly affected by processes not included in our 
models (shocks, X-rays). 

First, we have to calibrate qi. For that, we have pro¬ 
ceeded as follows: we have built the function qi(T e g , logg) 
using the same 2D interpolation method as for T e g (see 
Sect. 14.2.11) and then we have used the effective tem¬ 
peratures and gravities (relations T e g-ST and logg-ST) 
for each luminosity class to derive the ionising fluxes of 
dwarfs, giants and supergiants. Once again, we have used 
both theoretical and observational T e g scales to investi¬ 
gate the effects of such relations on the final values of qi 
(see Tables d to inj. Once obtained, these relations have 
been used with the relations R-ST to derive the calibra¬ 
tion Qi-spectral type. The results using the theoretical 
T e g scale are displayed in Fig. El] Figure ITU and IT1 show 
the theoretical and observational Qo’s and Qi’s and the 
differences between them. Although for Lyman fluxes the 
differences remain relatively small (lower than 0.3 dex), 
using different T e g scales can have a deep impact on the 
He I ionising fluxes, especially for late type supergiants for 
which Q\ may be changed by a factor 6 ! 

5.3.2. Comparisons with previous predictions 

Due to the presence of metals in the model atmospheres, 
ionising fluxes of O stars are modified: new bound-free 
opacities increase the blocking of flux at short wavelengths 
and the numerous metallic lines change the very shape of 
the spectral energy distribution. 

Giveon et al. ( 2002b ) and Morisset et al. (1200411 have 
shown that the new generation of atmosphere models in- 
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Spectral Type 

Fig. 13. Ionising fluxes as a function of spectral type for 
dwarfs (triangles), giants (squares) and supergiants (cir¬ 
cles) derived using the theoretical T e g scale. 

chiding a non-LTE approach, winds and line-blanketing 
produced more accurate ionising fluxes compared to other 
previous less sophisticated models. They used the new 
SEDs to compute photoionisation models and to build ex¬ 
citation diagrams (defined by the ratio of mid-IR nebular 
lines of Ne, Ar and S). The direct comparison of such dia¬ 
grams with observations of compact H II regions with ISO 
showed the improvement of the new generation models. 

Smith et al. computed new ionisation fluxes for 

O stars using WM-BASIC (Pauldrach et al. 1200111. Lanz & 
Hubeny (1200211 built a grid of plane parallel line-blanketed 
non-LTE models with TLUSTY and provided H and Hel 
ionising fluxes. Figures ED and El show the comparison 
between our ionising fluxes as a function of T e ff for dwarfs 
and supergiants with those of Smith et al. cop and the 
TLUSTY grid. Concerning H ionising fluxes, one notes 
that CMFGEN and TLUSTY agree very well (within 
0.1 dex). WM-BASIC models also show similar fluxes for 
supergiants and for the hottest dwarfs, while logqo are 
slightly lower (up to 0.3 dex) for dwarfs with T e ff < 40000 
K. The different gravity between our CMFGEN dwarfs 
(logg ~ 3.9) and the WM-BASIC dwarfs (logg = 4.0) 
may be partly responsible for this difference. However, 
since the difference in log qo are equivalent to a shift in T e g 
by less than 2000 K (see Fig. EH- we may speculate that a 
different treatment of line-blanketing between CMFGEN 
and WM-BASIC may be the main reason for the discrep¬ 
ancy. Concerning the Hel ionising fluxes, CMFGEN and 
TLUSTY are also in good agreement, which may be sur¬ 
prising given that TLUSTY models do not include winds. 
However, the largest wind effects are seen in the Hell con¬ 



Fig. 14. Lyman continuum photon fluxes derived using 
the theoretical (bold lines) and observational (thin lines) 
T e g scales. Solid lines are for dwarfs, dashed lines for gi¬ 
ants and dot-dashed lines for supergiants. The upper panel 
shows the difference between theoretical and observational 
results. 

tinuum ('Hillier H987alll987bl Gabler et al. 119891 Schaerer 
& de Koter 11997 ). Again, the WM-BASIC Hel ionising 
fluxes are lower for cool dwarfs (~ 0.3 dex) and supergiants 
(~ 0.1 dex). This results in softer spectra for WM-BASIC 



Fig. 15. Same as Fig.Elbut for He I ionising fluxes. 
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Fig. 16. Comparison between CMFGEN, TLUSTY and 
WM-BASIC ionising fluxes for dwarfs. TLUSTY models 
are from Lanz & Hubeny 1200211 and WM-BASIC models 
from Smith et al. m . Note that the latter have been 
estimated from the original SEDs and not from the SEDs 
re-binned for the inclusion in Starburst99. The CMFGEN 
fluxes are taken from Tables ^ to ED 

models compared to both CMFGEN and TLUSTY mod¬ 
els for late dwarfs and supergiants. The exact reasons 
for this difference leading to an apparent relative over¬ 
estimate of line blanketing remain unclear but may be 
likely attributed to different methods used to treat line¬ 
blanketing (opacity sampling versus direct inclusion with 
super-level). The behaviour of several atmosphere mod¬ 
els including CMFGEN, TLUSTY and WM-BASIC have 
been discussed extensively by Morisset et al. (IWCTl . In 
short, the CMFGEN and WM-BASIC models do a rea¬ 
sonable job when compared to observed fine structure line 
ratios of Galactic H II regions. However, at present the 
nebular analysis does not allow to favour CMFGEN or 
WM-BASIC models. New tailored studies of individual 
nebulae (following e.g. the steps of Oey et al. l?0DDt using 
these fully blanketed non-LTE models would be required 
to provide more strongest tests. 

The relations qi(T e g ) allow a direct comparison of at¬ 
mosphere models, but the quantities of interest for a num¬ 
ber of studies (H II regions, nebular analysis...) are usually 
the integrated ionising fluxes Q, calibrated as function of 
spectral type. Such relations ( Qi~ST ) depend not only 
on the model ingredients, but also on the effective tem¬ 
perature scale and on the R-ST relation. Figures IB IB 
and m show how our ionising fluxes - spectral type rela¬ 
tions derived using both the theoretical and observational 
effective temperature scales compare to the previous cal- 


Fig. 17. Same as Fig.^Jfor supergiants. 


ibrations of Vacca et al. and Schaerer & de Koter 

from CoStar models). 

Depending on which T e g scale is used, two kinds of con¬ 
clusions can be drawn. First, the theoretical ionising fluxes 
are lower than in the previous calibrations. Compared to 
the results of Yacca et al. GMJ, our Qo’s are 0.20 to 
0.80 dex lower for dwarfs, 0.25 to 0.55 lower for giants 
and 0.30 to 0.55 dex lower for supergiants. For Q i, the 
reduction is 0.1 to 1.6(!) dex for LC V, 0.2 to 0.8 dex for 
LC III and 0.25 to 0.8 dex for LC I. Such reductions are 
due to two distinct effects: first, the better treatment of 
line-blanketing compared to the CoStar models and the 
inclusion of both metals and winds compared to the re¬ 
sults of Vacca et al. cssa, which modifies the relations 
q i(T’eff); second, new relations T e ff-ST and R ST are used 
to compute Q,;. These effects render the CMFGEN spec¬ 
tra softer than the CoStar spectra (logQi/Qo lower by ~ 
0.2 to 1.0 dex). The Vacca et al. (1199611 predictions are 
usually harder than the CMFGEN SEDs, except for the 
earliest spectral types. 

Now if we use the T e g scales derived from the results of 
detailed spectroscopic analysis of individual massive stars, 
the ionising fluxes of stars with spectral types earlier than 
06 are basically unchanged compared to the results based 
on the “theoretical” T e g scale. For later spectral types the 
reduction of Qi for a given ST is much weaker, Qq and 
Q i being essentially similar to the Costar and Vacca et 
al. (119961) results for late supergiants. Quantitatively, the 
reduction of Qo is 0.25 to 0.50 dex (0.20 to 0.35 dex) for 
dwarfs (giants and supergiants). Q± is much less reduced 
with the observational scale (0.15 to 1.0 dex for dwarfs 
and 0.15 to 0.45 dex for giants) and are in fact similar 
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Fig. 18. Comparison between the present ionising fluxes 
of dwarfs and the results of Vacca et al. 119961 and Schaerer 
& de Koter 0223 CoStar models). 

—within 0.23 dex— to the Vacca et al. cn values for 
late supergiants. 

What do we learn from this? The obvious conclusion 
is that the choice of the underlying effective tempera¬ 
ture scale is important to determine the ionising fluxes 
for a given spectral type and luminosity class. As men¬ 
tioned in Sect. IP both the theoretical and observa¬ 
tional scales have their advantages and drawbacks. At 
present, we cannot say that one is more relevant than the 
other, which leaves us with an uncertainty on the ionising 
fluxes. However, both types of T e g - scales are significantly 
cooler than the ones based on plane-parallel H He models. 
Besides this, it is also worth noting that ionising fluxes 
are very sensitive to T e g (especially at low values) so that 
the difference we find between calibrations of QiS using 
different T e g - scales also puts forward the fact that the 
natural spread in T e g for a given spectral type (see Figs. 
El to HU leads to an uncertainty on QiS. And the error on 
the spectral type (of the order 0.5 unit) also translates to 
an error on the ionising fluxes (which amounts to a factor 
of ~ 2 on Q i for late spectral types). Thus, the present 
results are an improvement over previous analysis and cal¬ 
ibrations, but it should be reminded that they still suffer 
from uncertainties. 

6. Discussion: metallicity effects 

The present study is restricted to the case of a solar metal¬ 
licity for atmosphere models. First, we want to highlight 
that we have used a solar composition for all luminosity 
classes. This assumption is certainly good for dwarfs, but 



ST 

Fig. 19. Same as Fig. for giants. 



ST 

Fig. 20. Same as Fig.^Jfor supergiants. 

may be too crude for supergiants which usually show hints 
of He enrichment and modified CNO abundances (e.g. 
Herrero et al. Walborn et al. mi- These changes 

are likely different from star to star since they depend on 
rotational velocities (Meynet & Maeder l2000ll . Hence, non 
solar He abundance may slightly affect the strength of the 
classification lines and introduce a dispersion in ST for a 
given T e ff. Also, although Iron dominates the blanketing 
effect in terms of lines, CNO elements are important for 


log Q, log Qj 
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Table 1. Stellar parameters as a function of spectral types for luminosity class V stars obtained with the theoretical 
T e ff scales. 


ST 

Teff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 

log T© 

R 

[R©] 

M-spec 

[M©] 

log g 0 

[cm -2 s -1 ] 

log qi 

[cm" 2 s” 1 ] 

log Qo 
[s- 1 ] 

log Qi 

[s- 1 ] 

3 

44616 

3.92 

-5.79 

-4.03 

5.83 

13.84 

58.34 

24.56 

23.96 

49.63 

49.02 

4 

43419 

3.92 

-5.50 

-3.95 

5.68 

12.31 

46.16 

24.50 

23.86 

49.47 

48.83 

5 

41540 

3.92 

-5.21 

-3.82 

5.51 

11.08 

37.28 

24.38 

23.71 

49.26 

48.59 

5.5 

40062 

3.92 

-5.07 

-3.71 

5.41 

10.61 

34.17 

24.28 

23.59 

49.11 

48.42 

6 

38151 

3.92 

-4.92 

-3.57 

5.30 

10.23 

31.73 

24.15 

23.39 

48.96 

48.19 

6.5 

36826 

3.92 

-4.77 

-3.47 

5.20 

9.79 

29.02 

24.04 

23.17 

48.80 

47.93 

7 

35531 

3.92 

-4.63 

-3.36 

5.10 

9.37 

26.52 

23.91 

22.89 

48.63 

47.62 

7.5 

34419 

3.92 

-4.48 

-3.27 

5.00 

8.94 

24.15 

23.76 

22.49 

48.44 

47.18 

8 

33383 

3.92 

-4.34 

-3.18 

4.90 

8.52 

21.95 

23.65 

22.08 

48.29 

46.73 

8.5 

32522 

3.92 

-4.19 

-3.10 

4.82 

8.11 

19.82 

23.50 

21.58 

48.10 

46.18 

9 

31524 

3.92 

-4.05 

-3.01 

4.72 

7.73 

18.03 

23.34 

21.21 

47.90 

45.77 

9.5 

30488 

3.92 

-3.90 

-2.91 

4.62 

7.39 

16.46 

23.04 

20.59 

47.56 

45.12 


Table 2. Same as Tabled for luminosity class III stars. 


ST 

T e ff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 

log 4 

R 

[R©] 

M-spec 

[M©] 

log g 0 

[cm- 2 s” 1 ] 

log gi 

[cm" 2 s” 1 ] 

log Qo 
[s” 1 ] 

log Qi 
[s- 1 ] 

3 

42942 

3.77 

-6.13 

-3.92 

5.92 

16.57 

58.62 

24.49 

23.82 

49.71 

49.04 

4 

41486 

3.73 

-5.98 

-3.82 

5.82 

15.83 

48.80 

24.41 

23.71 

49.59 

48.89 

5 

39507 

3.69 

-5.84 

-3.67 

5.70 

15.26 

41.48 

24.29 

23.53 

49.44 

48.68 

5.5 

38003 

3.67 

-5.76 

-3.56 

5.63 

15.13 

38.92 

24.19 

23.39 

49.33 

48.53 

6 

36673 

3.65 

-5.69 

-3.45 

5.56 

14.97 

36.38 

24.08 

23.23 

49.22 

48.37 

6.5 

35644 

3.63 

-5.62 

-3.37 

5.49 

14.74 

33.68 

23.99 

23.07 

49.11 

48.19 

7 

34638 

3.61 

-5.54 

-3.28 

5.43 

14.51 

31.17 

23.89 

22.85 

49.00 

47.96 

7.5 

33487 

3.59 

-5.47 

-3.18 

5.36 

14.34 

29.06 

23.79 

22.55 

48.89 

47.64 

8 

32573 

3.57 

-5.40 

-3.10 

5.30 

14.11 

26.89 

23.68 

22.21 

48.76 

47.29 

8.5 

31689 

3.55 

-5.32 

-3.02 

5.24 

13.88 

24.84 

23.57 

21.90 

48.64 

46.97 

9 

30737 

3.53 

-5.25 

-2.93 

5.17 

13.69 

23.07 

23.37 

21.43 

48.42 

46.49 

9.5 

30231 

3.51 

-5.18 

-2.88 

5.12 

13.37 

21.04 

23.27 

21.18 

48.30 

46.22 


Table 3. Same as Tabled for luminosity class I stars. 


ST 

T e ff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 

log 4 

R 

[R©] 

M-spec 

[Mo] 

log g 0 

[cm -2 s" 1 ] 

log gi 

[cm" 2 s -1 ] 

log Qo 
[s- 1 ] 

log Qi 

K 1 } 

3 

42551 

3.73 

-6.35 

-3.89 

6.00 

18.47 

66.89 

24.47 

23.79 

49.79 

49.11 

4 

40702 

3.65 

-6.34 

-3.76 

5.94 

18.91 

58.03 

24.38 

23.64 

49.71 

48.98 

5 

38520 

3.57 

-6.33 

-3.60 

5.87 

19.48 

50.87 

24.25 

23.45 

49.61 

48.81 

5.5 

37070 

3.52 

-6.33 

-3.48 

5.82 

19.92 

48.29 

24.15 

23.31 

49.53 

48.69 

6 

35747 

3.48 

-6.32 

-3.38 

5.78 

20.33 

45.78 

24.04 

23.14 

49.44 

48.54 

6.5 

34654 

3.44 

-6.31 

-3.29 

5.74 

20.68 

43.10 

23.95 

22.97 

49.36 

48.39 

7 

33326 

3.40 

-6.31 

-3.17 

5.69 

21.14 

40.91 

23.84 

22.69 

49.27 

48.12 

7.5 

31913 

3.36 

-6.30 

-3.04 

5.64 

21.69 

39.17 

23.68 

22.32 

49.14 

47.77 

8 

31009 

3.32 

-6.30 

-2.96 

5.60 

22.03 

36.77 

23.56 

21.97 

49.03 

47.45 

8.5 

30504 

3.28 

-6.29 

-2.91 

5.58 

22.20 

33.90 

23.48 

21.75 

48.95 

47.22 

9 

29569 

3.23 

-6.29 

-2.82 

5.54 

22.60 

31.95 

23.31 

21.22 

48.80 

46.71 

9.5 

28430 

3.19 

-6.28 

-2.70 

5.49 

23.11 

30.41 

23.17 

21.05 

48.69 

46.57 


the blocking of radiation through their bound-free tran¬ 
sitions usually close to the He II edge so that, although 
the sum of C, N and O abundances are constant during 
CNO processing, a change in the abundance of each el¬ 
ement may modify the global blanketing effect. Such an 
effect should be further studied by new analyses of su¬ 


pergiants. Second, the present work does not address the 
important question of the dependence of the fundamental 
parameters of O stars on global metallicity. Nonetheless, 
several indications of the general trend exist. 

From the modelling side, Kudritzki 2()()~2; and Mokiem 
et al. mm have investigated the effect of a change of 
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Table 4. Stellar parameters as a function of spectral types for luminosity class V stars obtained with the observational 
T e ff scales. 


ST 

Teff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 

lo s4 

R 

[Re] 

M-spec 

[M©] 

log g 0 

[cm -2 s 1 ] 

log qi 

[cm" 2 s” 1 ] 

log Qo 
[s- 1 ] 

log Qi 

[s- 1 ] 

3 

44852 

3.92 

-5.79 

-4.05 

5.84 

13.80 

57.95 

24.57 

23.98 

49.64 

49.04 

4 

42857 

3.92 

-5.50 

-3.91 

5.67 

12.42 

46.94 

24.47 

23.82 

49.44 

48.79 

5 

40862 

3.92 

-5.21 

-3.77 

5.49 

11.20 

38.08 

24.34 

23.66 

49.22 

48.54 

5.5 

39865 

3.92 

-5.07 

-3.70 

5.41 

10.64 

34.39 

24.27 

23.57 

49.10 

48.41 

6 

38867 

3.92 

-4.92 

-3.62 

5.32 

10.11 

30.98 

24.20 

23.46 

48.99 

48.26 

6.5 

37870 

3.92 

-4.77 

-3.55 

5.23 

9.61 

28.00 

24.13 

23.34 

48.88 

48.09 

7 

36872 

3.92 

-4.63 

-3.47 

5.14 

9.15 

25.29 

24.04 

23.17 

48.75 

47.88 

7.5 

35874 

3.92 

-4.48 

-3.39 

5.05 

8.70 

22.90 

23.94 

22.98 

48.61 

47.64 

8 

34877 

3.92 

-4.34 

-3.30 

4.96 

8.29 

20.76 

23.82 

22.68 

48.44 

47.30 

8.5 

33879 

3.92 

-4.19 

-3.22 

4.86 

7.90 

18.80 

23.69 

22.23 

48.27 

46.81 

9 

32882 

3.92 

-4.05 

-3.13 

4.77 

7.53 

17.08 

23.52 

21.69 

48.06 

46.22 

9.5 

31884 

3.92 

-3.90 

-3.04 

4.68 

7.18 

15.55 

23.39 

21.31 

47.88 

45.80 


Table 5. Same as Table 0 for luminosity class III stars. 


ST 

T e ff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 


R 

F®] 

M-spec 

[M©] 

log g 0 

[cm- 2 s” 1 ] 

log gi 

[cm" 2 s” 1 ] 

log Qo 
[s” 1 ] 

log Qi 
[s- 1 ] 

3 

44537 

3.77 

-6.13 

-4.03 

5.96 

16.19 

55.95 

24.57 

23.94 

49.77 

49.14 

4 

42422 

3.73 

-5.98 

-3.88 

5.85 

15.61 

47.43 

24.47 

23.78 

49.64 

48.95 

5 

40307 

3.69 

-5.84 

-3.73 

5.73 

15.07 

40.43 

24.34 

23.61 

49.48 

48.75 

5.5 

39249 

3.67 

-5.76 

-3.65 

5.67 

14.82 

37.35 

24.28 

23.51 

49.40 

48.63 

6 

38192 

3.65 

-5.69 

-3.57 

5.61 

14.58 

34.53 

24.21 

23.41 

49.32 

48.52 

6.5 

37134 

3.63 

-5.62 

-3.49 

5.54 

14.36 

31.96 

24.13 

23.28 

49.23 

48.38 

7 

36077 

3.61 

-5.54 

-3.40 

5.48 

14.14 

29.59 

24.04 

23.16 

49.13 

48.24 

7.5 

35019 

3.59 

-5.47 

-3.32 

5.42 

13.93 

27.45 

23.94 

22.96 

49.01 

48.03 

8 

33961 

3.57 

-5.40 

-3.23 

5.35 

13.74 

25.49 

23.82 

22.68 

48.88 

47.74 

8.5 

32904 

3.55 

-5.32 

-3.13 

5.28 

13.55 

23.68 

23.70 

22.32 

48.75 

47.37 

9 

31846 

3.53 

-5.25 

-3.04 

5.21 

13.38 

22.04 

23.61 

22.02 

48.65 

47.06 

9.5 

30789 

3.51 

-5.18 

-2.94 

5.15 

13.22 

20.55 

23.39 

21.50 

48.42 

46.53 


Table 6. Same as Table 0 for luminosity class I stars. 


ST 

T e ff 

[K] 

log Qspec 

[cm s -2 ] 

M v 

BC 


R 

[Re] 

M-spec 

[Mo] 

log g 0 

[cm" 2 s" 1 ] 

log gi 

[cm" 2 s -1 ] 

log Qo 
[s- 1 ] 

log Qi 

K 1 } 

3 

42233 

3.73 

-6.35 

-3.87 

5.99 

18.56 

67.53 

24.46 

23.77 

49.78 

49.09 

4 

40422 

3.65 

-6.34 

-3.74 

5.93 

18.99 

58.54 

24.36 

23.62 

49.70 

48.96 

5 

38612 

3.57 

-6.33 

-3.61 

5.87 

19.45 

50.72 

24.25 

23.46 

49.62 

48.82 

5.5 

37706 

3.52 

-6.33 

-3.54 

5.84 

19.70 

47.25 

24.20 

23.38 

49.58 

48.76 

6 

36801 

3.48 

-6.32 

-3.46 

5.81 

19.95 

44.10 

24.14 

23.29 

49.52 

48.67 

6.5 

35895 

3.44 

-6.31 

-3.39 

5.78 

20.22 

41.20 

24.07 

23.18 

49.46 

48.58 

7 

34990 

3.40 

-6.31 

-3.31 

5.75 

20.49 

38.44 

24.00 

23.06 

49.41 

48.47 

7.5 

34084 

3.36 

-6.30 

-3.24 

5.72 

20.79 

36.00 

23.89 

22.89 

49.31 

48.31 

8 

33179 

3.32 

-6.30 

-3.16 

5.68 

21.10 

33.72 

23.81 

22.70 

49.25 

48.14 

8.5 

32274 

3.28 

-6.29 

-3.08 

5.65 

21.41 

31.54 

23.74 

22.56 

49.19 

48.00 

9 

31368 

3.23 

-6.29 

-2.99 

5.61 

21.76 

29.63 

23.65 

22.21 

49.11 

47.67 

9.5 

30463 

3.19x 

-6.28 

-2.91 

5.57 

22.11 

27.83 

23.52 

21.87 

49.00 

47.35 


the metal content on the spectral energy distribution of O 
dwarfs using CMFGEN models. They found that H ion¬ 
ising fluxes are essentially unchanged when Z is varied 
between twice and one tenth the solar content. They ar¬ 
gue that the redistribution of the flux blocked by metals 
at short wavelengths takes place within the Lyman con¬ 


tinuum, which explains the observed behaviour. However, 
they show that the SEDs are strongly modified below ~ 
450 A, spectra being softer at higher metallicity (see also 
Sect. l5.fl.Tl) . Morisset et al. 1200411 have computed various 
WM-BASIC models at different metallicities and showed 
how Z affected the strength of mid-IR nebular lines emit- 
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ted in compact H II regions. The softening of the SEDs 
when metallicity increases is crucial to understand the be¬ 
haviour of observed excitations sequences. 

Concerning the effect of metallicity on T e g scales, 
Mokiem et al. 1200411 have shown that for a given T e g , 
spectral types vary within one subclass when Z is de¬ 
creased from 2 to 0.1 Zq. This boils down to a higher ef¬ 
fective temperature at low metallicity than at solar metal¬ 
licity. We reach the same conclusion from the study of sev¬ 
eral test models for dwarfs with Z = 1/8 Zq (see Martins 
l2004btl . We estimate that the reduction of the T e g scale 
(compared to pure H He results) at this metallicity is 
roughly half the reduction obtained at solar metallicity 
(see also Martins et al. rna . 

Observationally, Massey et al. 1200411 derived effective 
temperatures of a sample of O stars in the Magellanic 
Clouds by means of models computed with the code 
FASTWIND (Santolaya-Rey et al. They found 

lower T e g than Vacca et al. (THTO but higher than those of 
Galactic O stars, in good agreement with our results. They 
estimate that effective temperatures of early to mid O 
type MC objects are 3000 to 4000 K hotter than Galactic 
counterparts. Previously, Crowther et al. 12002 |l computed 
CMFGEN models and showed that extreme O supergiants 
in the Magellanic Clouds were cooler by 5000 to 7500 K 
compared to the Vacca et al. dsns calibration. Bouret 
et al. mm also confirmed the reduction of effective 
temperatures in SMC O dwarfs using CMFGEN models. 
However, they derived temperatures in agreement our T e g 
scale for Galactic stars, which may be surprising given the 
lower metallicity of the SMC. Even more recently, Heap et 
al. iffiTO submitted a paper in which they derive the stel¬ 
lar properties of a large sample of SMC stars. Again, the 
effective temperatures are not inconsistent with those of 
Galactic stars. They also provide bolometric corrections 
and ionising fluxes which are in better agreement with 
our new calibration than with the older pure H He plane- 
parallel ones. Given the dispersion in the data points, a 
conclusion as regards the metallicity dependence is not 
possible. 

The above discussion shows that although hints on the 
metallicity dependence of stellar parameters of O stars 
exists, a quantitative and systematic study remains to be 
carried out. 


7. Concluding remarks 

We have presented new calibrations of stellar parameters 
of O stars as a function of spectral type based on atmo¬ 
sphere models computed with the code CMFGEN (Hillier 
& Miller IlOQTi ). We have build a grid of models spanning 
the range of spectral types and luminosity classes of O 
stars from which calibrations have been derived through 
2D interpolations. The main improvement of such rela¬ 
tions over previous ones is the inclusion of line-blanketing 
and winds in the non-LTE atmosphere models. Our main 
results are the following: 


o We have derived two types of effective temperature 
scales: a theoretical one based uniquely on our grid of 
models (approach similar to paper I) and an observa¬ 
tional one from the results of detailed spectroscopic 
analysis of individual O stars with line blanketed non- 
LTE models including winds. We confirm the now well 
established fact that line-blanketing leads to cooler ef¬ 
fective temperature scales compared to the widely used 
relations of Vacca et al. n™ based on plane-parallel 
pure H He models. The adoption of a better photo- 
spheric structure in our models leads to slightly cooler 
(~ 500 K) theoretical effective temperatures compared 
to our first results (Martins et al. 2002 :. Our theoret¬ 
ical T e ff-ST relations are cooler by 2000 to 8000 K 
being larger for earlier spectral types and lower lu¬ 
minosity classes — compared to Vacca et al. mu 
calibrations. The theoretical T e g scales are consistent 
with observational relations for early type dwarfs and 
supergiants and are slightly cooler for late type spec¬ 
tral types (by up to 2000 K for supergiants). This may 
be due to too high mass loss rates in our models. For 
giants, the theoretical relation is systematically cooler 
by ~ 1000 K. We estimate the uncertainty of our the¬ 
oretical effective temperatures to be between 1000 and 
2000 K for a given spectral type. This is due to the in¬ 
trinsic dispersion of temperatures for O stars of a given 
spectral type and to the uncertainties of the models 
(hydrodynamic structure, adopted parameters), 
o Bolometric corrections are reduced compared to the 
calibrations of Vacca et al. dmu, the reduction be¬ 
ing larger for supergiants and for the earliest spectral 
types. BC’s estimated using the theoretical T e g scales 
are 0.40 to 0.50 (resp. 0.35 to 0.60, 0.40 to 0.65) mag 
lower for dwarfs (resp. giants, supergiants), while those 
estimated using the observational relations are 0.30 to 
0.50 (resp. 0.30 to 0.50, 0.20 to 0.65) mag lower for 
luminosity class V (resp. Ill, I). 
o Luminosities are reduced by 0.20 to 0.35 dex for 
dwarfs, 0.25 dex for giants and 0.25 to 0.35 dex for 
supergiants compared to the calibrations of Vacca et 
al. Ill 996|) with little difference (< 0.1 dex) between re¬ 
sults using theoretical and observational relations. The 
reduction of luminosity is independent of spectral type 
for giants and supergiants and is slightly larger for late 
type than for early type dwarfs, 
o For a given spectral type, ionising fluxes are reduced 
due to the effect of line-blanketing on the SED (block¬ 
ing of flux). Compared to the Vacca et al. 11199611 val¬ 
ues and using our theoretical T e g scale, we find inte¬ 
grated Lyman ionising fluxes (Qo) 0.20 to 0.80 dex 
lower for dwarfs, the larger difference being at late 
spectral types. The reduction is 0.25 to 0.55 dex (resp. 
0.30 to 0.55 dex) for giants (resp. supergiants). If we 
use the observational effective temperature scale, the 
reduction are 0.25 to 0.50 dex for dwarfs and 0.20 to 
0.35 dex for giants and supergiants. He I ionising fluxes 
are also reduced, the difference between theoretical and 
observational results being larger. 
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For a given T e g, our g*’s agree well with the TLUSTY 
grid OSTAR2002 (Lanz & Hubenv [20021) . but for late 
spectral types, they are larger by up to 0.3 dex com¬ 
pared to the results of Smith et al. 1200211 based on 
WM-BASIC models. This shows that current atmo¬ 
sphere codes are not fully consistent and indicates that 
the typical uncertainty on current ionising fluxes per 
unit area may still be up to a factor of ~ 2 for a given 
effective temperature. 

Our results should be tested further when more analy¬ 
sis of individual stars are available and cover the whole 
range of spectral types and luminosity classes. Despite 
some remaining uncertainties on the T e g scale, our results 
should represent a significant improvement over previous 
calibrations, given the detailed treatment of non-LTE line¬ 
blanketing in the expanding atmospheres of massive stars. 
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